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The rates oi + n ^ e~ + p and + p ^ + n are important for understanding the dynamics 
of supernova explosion and the production of heavy elements in the supernova environment above 
the protoneutron star. Observations and theoretical considerations suggest that some protoneutron 
stars may be born with strong magnetic fields. In a previous paper we calculated the above rates 
in supernova environments with magnetic fields up to ~ 10^'' G assuming that the nucleon mass 
miv is infinite. We also applied these rates to discuss the implications of such strong fields for 
supernova dynamics. In the present paper we take into account the effects of a finite niM and 
develop a numerical method to recalculate the above rates in similar environments. This method is 
accurate to C(l/mjv) and such an accuracy is required for application to supernova nucleosynthesis. 
We show that our results have the correct behavior in the limit of high neutrino energy or small 
magnetic field. Based on comparison of our results with various approximations, we recommend 
efficient estimates of the above rates for use in models of supernova nucleosynthesis in the presence 
of strong magnetic fields. 

PACS numbers: 25.30.Pt, 26.50.+X, 97.60.Bw 



I. INTRODUCTION 

The processes 

i^e + n ^ e~ + p, (la) 
De+P ^ e+ + n, (lb) 

play important roles in supernovae. A supernova is initiated by the collapse of a stellar core, which leads to the 
formation of a protoneutron star. Nearly all the gravitational binding energy of the protoneutron star is emitted in 
i/fi, j^e, Vfi, Vr, and Vr, some of which would interact to heat the material above the protoneutron star. The forward 
neutrino absorption processes in Eq. provide the dominant heating mechanism, which is counteracted by cooling 
of the material through the reverse n7utrino emission processes. In a prevalent paradigm Q , supernova explosion 
is determined by the competition between these heating and cooling processes. These processes also interconvert 
neutrons and protons, thereby setting the neutron-to-proton ratio of the material above the protoneutron star 
This ratio is a key parameter that governs the production of heavy elements during the ejection of this material 
[^Q. Thus, accurate rates of the processes in Eq. ^ arc important for understanding supernova dynamics and 
nucleosynthesis. 

Observations and theoretical considerations indicate that protoneutron stars with magnetic fields of ^ 10^^ G may 
be formed. The rates of the abo ve p rocesses in such strong fields have been studied in the literature with various 
approximations [slIR H. lol ITol UlTIT^. In our previous work |0, we used the Landau wave functions of and 
derived a set of simple and consistent formulas to calculate the rates of the processes in Eq. (^Q in the presence of 
strong magnetic fields. We also applied these rates to discuss the implications of such fields for supernova dynamics. 
However, all the calculations in the literature, including our previous work, were mostly carried out to 0(1), the 
zeroth order in I/toat with mjv being the nucleon mass. None of them included both the effects of nucleon recoil 
and weak magnetism, which are of ©(I/ttiat) and known to be important for the conditions in supernovae 14J. For 
modeling the production of heavy elements during the ejection of the material from the protoneutron star, an accuracy 
of ~ 1% for the rates of the processes in Eq. iQJ is required to determine precisely the neutron-to-proton ratio in 
the material P|. To achieve such an accuracy, the ©(l/m^v) effects on these rates must be taken into account. In 
this paper we recalculate these rates using the respective Landau wave functions of and protons and including the 
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©(l/mjv) corrections from both nucleon recoil and weak magnetism. Our goal is to identify the important factors in 
computing accurate rates of the processes in Eq. fQ) for application to supernova nucleosynthesis in the presence of 
strong magnetic fields. 

This paper is organized as follows. The energies and the wave functions of the relevant particles in magnetic fields 
are discussed in Sec.m The cross sections of the neutrino absorption processes in Eq. Q and the differential reaction 
rates of the reverse neutrino emission processes are derived to 0{l/mN) in Sec. IIIII The rates of these processes in 
supernova environments with magnetic fields of ^ 10^^ G are calculated and discussed in Sec. IIVI Conclusions are 
given in Sec. El 



II. PARTICLE ENERGIES AND WAVE FUNCTIONS IN MAGNETIC FIELDS 



The importance of magnetic field effects can be gauged from the energy scale 

^ Xl/2 



Vei3 = 7.69(^^^j MeV, (2) 

where e is the charge of e"*" and B is the field strength. However, there is no detailed knowledge of magnetic fields 
in supernovae. Observations indicate that neutron stars may have B 10^^ G long after their birth in supernovae 
[lalla . ll7j . This suggests that at least B ^ 10^^ G can be generated during the formation of some protoneutron stars. 
A recent theoretical model suggests that B ^ 10^^ G may be produced near the surface of a protoneutron star Ts"!. 
An upper limit of B ^ 10^^ G can be estimated for such a star by equating the magnetic energy to its gravitational 
binding energy . To explore the effects of strong magnetic fields on the rates of the processes in Eq. ^ , we consider 
that B ~ 10^^ G may exist in the region of interest to supernova nucleosynthesis, which lies well below 10^ cm from 
the protoneutron star 0. For such fields, the associated energy scale is much smaller than the mass of the W boson 
Mw = 80 GeV. So there will be no change in the description of the weak interaction that is involved in the processes 
in Eq. (Q. On the other hand, the energy scale for B ^ 10^^ G is larger than the temperature (T ~ f MeV) of the 
material above the protoneutron star and comparable to the typical neutrino energy (E^, ~ 10 MeV). Thus, magnetic 
field effects on energy levels of charged particles (e^ and p) will be important. Furthermore, B ^ 10^^ G will induce 
polarization of nucleon spin at the level of eB/rriNT 10~^. This is significant due to parity violation of weak 
interaction and should also be taken into account. 

We discuss the energy levels and the corresponding wave functions of all the relevant particles in this section. We 
assume a uniform magnetic field B in the positive z-direction, for which the vector potential is 

A=i-^By,^Bx,0). (3) 
All the wave functions will be given in Dirac-Pauli representation. 



A. Electron and positron 

The motion of along the z-axis is not affected by the magnetic field, but the motion in the xy-plane is quantized 
into Landau levels with energies (see, e.g., Ref. [2^) 

Ee = ^/ml + kj, + 2neeB, (4) 

where nie is the rest mass of e^, k^z is the z-component of the momentum, and Ue is an integer quantum number (i.e., 
Ue — 0, 1, 2, . . . ). For the in the initial states of the neutrino emission processes in Eq. Q), the relevant Ee is of 
the order of the temperature T ^ 1 MeV for the material above the protoneutron star, ft can be seen from Eqs. (O 
and I0J that these predominantly occupy the ground Landau level (rie = 0) for B ^ 10^^ G. In comparison, the 
in the final states of the neutrino absorption processes typically have E^ of the order of the neutrino energy Ei, ^ 10 
MeV. These can occupy excited Landau levels {rig > 1). 
The wave function of in cylindrical coordinates (^, 0, z) is 



(5) 
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where Se = 1 and —1 for spin up and down, respectively, rg is the quantum number labeUng the center of gyromotion 
in the xy-plane, and L is the linear size of the normalization volume. In Eq. Q, the spinor {U^-)s^ is 



{Ue-)s 





fce, e-"^/„ ,_i.,,(egeV2) 



(6) 



and 



\ 



I 

(me + £ ;e)/„,..,(e5gV2) 
-i^2n^eBe-"i'In,-\,r,_ {eB^"^ /2) 



(7) 



The special function /„,r(C) in the above equations is defined in Ref. |[2l| and can be calculated using the method 
given in Appendix IXI 

The wave function of e"*" is 



^2nL/eB 



-(f^e+)s., 



(8) 



where 



y/2Ee{E^+me) 



( W2neeBe-"t'In,-i^r^ {eB^^ 12) \ 

-fce./„..re(e^eV2) 



\ (me+£;e)/n..r.(eS^V2) / 



(9) 



and 



e+ Jse=-1 



^2E^{E^ + vne) 



\ 



-fce .e-^^/„ ,-i,,,(eBeV2) \ 
W2neeBIn^^r, {eB^^/2) 
-(me + £;e)e-'*i„,_i,,,(eBeV2) 
/ 



(10) 



Clearly, E,, does not depend on the quantum number in the wave functions. This leads to a degeneracy factor 

eBL^ 



El 



27r 



(11) 



for each Landau level of (see, e.g., Ref. ^3). Each level is further degenerate with respect to spin except for the 
ground level [{ipe-)s^=i = ('0e+)sc=-i = for "^e = 0]. This introduces an additional spin degeneracy factor (7„^, 
which is 1 for rig = and 2 for nf,>l. 



B. Proton 



Protons are nonrelativistic in the supernova environment of interest. For nonrelativistic with the same charge 
and spin as protons, expansion of Eq. Q to ©(l/mg) gives 

fc^, UpcB 

Ee.NR = me + + . 12 

The above equation already accounts for the contribution from the e+ magnetic moment of e/2me- Unlike e"*", protons 
have an anomalous magnetic moment of flp — 1.79 fiN in addition to the value (Xn = e/2mp expected for a spin-1/2 
point particle of charge e and mass rup. Taking this into consideration, we obtain the energies of the proton Landau 
levels as 

fcpz , ripcB 

Ep^mp + - 1 SpfipB, (13) 
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where symbols have similar meanings to those for e . For B ~ 10 G, eB /nip flpB '--^ 60 keV. The protons in the 



initial states of the absorption and emission processes in Eq. have Ep 



T ~ 1 MeV, and therefore, 



can occupy many Landau levels. By the correspondence principle, the quantum effects of the magnetic field on these 
protons are insignificant. However, the protons in the final states of the absorption and emission processes are 
less energetic, with typical recoil energies of Ep — uip ~ E'^/nip ^ 100 keV and ~ T'^ /nip ~ 1 keV, respectively. Thus, 
proper treatment of proton Landau levels is especially important for these processes. 
The proton wave function can be written as 



^2-KL/eB 



(14) 



where 



In,,r,{eBe /'2) 


(fcp^/2TOp)/„p,rp(eB$V2) 
\ -ie'*(yi^^/2mp)/„,-i.r,(eBeV2) / 



(15) 



and 



(f/p)s, = -l = 



^^/„^_i,,^(eSeV2) 



i{^2^^/2nip)I,,^^r,{eBe/'i) 
V-e^^(fcp./2mp)/„^_i,,^(eBeV2) / 



(16) 



Note that each proton Landau level is also degenerate with respect to the quantum number rp, but the spin degeneracy 
of the excited levels is lifted due to the contribution from the anomalous magnetic moment. 
The proton wave function contains terms of the form 



y/2TrL/eB 



-Ir.AeBe/2), 



which has the following properties: 



and 



The operator tt in the above equations is defined as 

TT = — iV — eA. 



(17) 



(18a) 



(18b) 



(19) 



Equations (|18a|l and Ijl8bp can be used to simplify the evaluation of the transition amplitudes for the processes in 
Eq. O 



C. Neutron 



Neutrons are also nonrelativistic in the supernova environment of interest, and their energy is 



En = "In 



2nir, 



where //„ = — 1.91//Ar is the neutron magnetic moment. The corresponding wave function to ©(I/toat) is 

i(k„-x-£;„t) 



(^n)s 



L3/2 



(20) 



(21) 
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where 



{Un)s 



( 



(fc„/2m„)cos9„ 
\ {kn/2mn) sine„e**" 



(22) 



and 



/ 



{Un) 



s„=-l — 



(fc„/2m„) sine„e 
V -(fc„/2mn)cose„ 



(23) 



In the above equations, 0„ and are the polar and azimuthal angles of the neutron momentum k„ in spherical 
coordinates. 



D. Neutrinos 



The neutrino energy is not affected by the magnetic field. For left-handed Ve with momentum kj^, the wave function 



IS 



where 



oi(k„-x-E„t) 



L3/2 



U„ = 



sin(e,/2) 
-cos(e^/2) 
-sin(e,/2) 

cos(ei./2) 



(24) 



(25) 



The azimuthal angle of k^, is taken to be = in the above equation. The wave function of right-handed Ve with 
the same momentum kj^ is 



-j(k^-x-£;„t) 



(26) 



where Ui, is the same as in Eq. ()25|l . 



III. CROSS SECTIONS AND DIFFERENTIAL REACTION RATES 



As discussed in Sec. i? ^ 10^^ G will not affect the weak interaction, which is still described by the effective 
four-fermion Lagrangian 



(27) 



where Gp — (292.8 GeV) ^ is the Fermi constant, dc is the Cabbibo angle (cos^ 0c = 0.95), the leptonic charged 



current L" is 



and the nucleonic current N" is 



L" = V;.7"(l-75)V'e 



(28) 



(29) 



In the above equations, 7", 75, and cr"'^ are the standard matrices describing fermionic transitions in Dirac-Pauli 
representation, and f = I, g = 1.26, and /2 = 3.7 are the nucleon form factors. [A more up-to-date value of g is 
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1.27 [23. This value is recommended for calculating the rates of the processes in Eq. for specific application to 
supernova nucleosynthesis.] The term involving /2 in Eq. H29|l represents weak magnetism and must be included for 

calculations to 0{l/mN)- The covariant derivative —iVp in this term preserves the gauge invariance and operates 
according to 



il)pO {^-i V f}^ ipn = [i~idp - eAp) V'p] Oipn - tppO (idpipn) ■, 



(30) 



where O is a constant matrix and Afj corresponds to the electromagnetic field {Aq = here). 

Based on the above description of the weak interaction, we derive below the cross sections of the neutrino absorption 
processes in Eq. ^ and the differential reaction rates of the reverse neutrino emission processes. We will include 
the magnetic field effects on particle energies and wave functions and focus on corrections of C(l/m7v) in both the 
transition amplitude and kinematics. Radiative corrections and the effect of the Coulomb field of the proton on the 
electron wave function are ignored for simplicity. (The Coulomb field will modify the Landau wave function of the 
electron, thus making the calculation much more complicated.) We propose an approximate treatment of these factors 
at the end of Sec. IrT^ 



A. Cross sections for neutrino absorption 



We first derive the cross section of i/g + n ^ e + p in detail. The transition matrix of this process is 



Gf cos 9c 
V2 



2m„ V / 



With the wave functions given in Sec. ^1 Eq- H31|l can be rewritten as 

Gf cos 9c eB 



- V2 27TL 
The amplitude 97ly^„ in Eq. is 



an. 



Cdc 
^0 

2m,; 



gj:wj_ •xj_ g-i(nc -rc -rip+rp) 



(Xp)t7°a"^C/„ - {kn)0Upa''^Un\ U,-^M - 75)1/.. } dc/) 



(31) 



(32) 



(33) 



where w = k„ + is the total momentum, the subscript _L denotes a vector in the xy-plane, and 

-1 



{Xp)i3 = 



oi(fep^z-£;pt)gi(rp-np)0 
y/2TTL/eB 



{idp - eAjs) ipp. 



(34) 



Evaluation of {Xp)p for /? = 1 and 2 (x and y) can be simplified by using Eqs. (|18a|l and Ijl8b|l [2^. 

The (5-functions in Eq. H32|l enforce conservation of energy and of momentum in the 2;-direction, for which both 
the neutron and the proton momenta must be taken into account in calculations to 0(l/mAr). For + n — s- +p 
occurring in the material above the protoneutron star, the neutron momentum is especially important as the typical 
value fc„ ~ -\/2m„r — 43(r/MeV)^/^ MeV is larger than the typical Ve momentum k^ — E^ ^ 10 MeV. To account 
for this, we average the cross section over the normalized thermal distribution function /„(k„,s„) for the neutrons 
and obtain 



^l!n''^=E / /n(kn,S„)d3fc„ ^ 



Sp , Tip .Vp 



Ldk 



'El \^ 



Ldkp 



1 \%. 



27r ^ / 27r L-^L-^ tL^ 



(35) 



where the superscript (1,-B) denotes the cross section to ©(I/toat) and in the presence of a magnetic field, r is the 
duration of the interaction, and 



In (k/i , S^) 



,-(£;„-m„)/T 



(36) 
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The summation and integration in Eq. I|35() must be treated as nested integrals. For example, the summation over s„ 
and the integration over k„ apply to not only /„(k„, s„) but also the subsequent terms that have implicit dependence 
on Sn and k„. The summation and integration in the equations below should be interpreted similarly. 
Using 

y 6{Ee + Ep~ En) 6{kez + ~ Kz - knz) d^K 



= m„ j d$„ (37) 

to integrate over the neutron momentum, we can rewrite Eq. (|35(l as 

i^ ^ g^ _ G% cos^ 9c mncB 1 



4tt (27rm„T)3/2 et'^B/T g-t^^B/T 



EE/ dfce. E E /dfcp. / e-^^-'-^/^W.^^d*™, 



(38) 



where 



is the reduced amplitude squared given explicitly in Appendix El It follows from Eq. 1)37(1 that £'„ and knz in the 
integrand in Eq. H38|l are determined in terms of other quantities by conservation of energy and of momentum in the 
^-direction. The integration region K oi J dkpz in Eq. (|38|l is also set by these conservation laws, which require 



E, + m„ + ^+ ^^^^+^J>^ - Sn^^nB = Vml + kl + 2n,eB 



'^'^-Sp-^pB. (40) 



The above equation can be rearranged into the form 



ak' + bkpz + c = > 0, (41) 



2m„ 



where 



2mpmn ' 



(42a) 



Thus, 



h = , (42b) 

c = y/ml + kl, + 2neeB - E^ - A - ^^"^ ^''"^ + - SpjlpB + s„/i„B, (42c) 

A = rUn — mp. (42d) 



^ ^ f(-oo,+cx)), iffe2<4ac, 
|(-oo,(fcp2)_]U[(fcp^)+,+oo), if62>4ac. 



where 



(fcp.)± 2^ • (44) 
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For !/e + p ^ e+ + n occurring in the material above the protoneutron star, the cross section can be written as 

(27r)3 ^ J 2n L-^L-^ tL^ 



'^ip ^ - X! X! f fpikpz,np,Sp)dkpz^ f ^ / _ ' "3' , (45) 



where 



and 



is the normahzed thermal distribution function for the protons. Using again the integration over the neutron momen- 
tum to get rid of the (5-functions, we obtain 

_ Gl cos2 9c rrin 1 - e-^^/™^^ 



E E / dfce. ^ J2 /dfcpJ e-(^''-"^)/^>V,,pd<I>„. (48) 
The reduced amplitude squared Wp^p in the above equation can be obtained from Wu^n by making the substitution 

The integration region K' of / dfcpz in Eq. H48|l is determined from energy and momentum conservation as in Eq. 140(1 
but with the above substitution implemented. 

For application to supernova neutrinos, it is useful to further average the cross sections in Eqs. ()38|) and (|48|) over 
the relevant normalized neutrino energy spectra fu{Ey) to obtain 



j (J,NME.)dE,, (50) 



where ai,N stands for aljn or cr^^p . A typical form of f^{E,j) adopted in the literature is 

1 E^ 

^•"^^"^ ^ TSF2Me^p[iE,/T^)-r,,] + V ^^^^ 
where and rj^, are constant parameters and 

FnM^ / 7 VXT'^^- (52) 

Jo exp (a; - 7^,, ) + 1 

For the neutrino energy spectra in Eq. 151() . the average neutrino energy is 

(E.) = ^,T^- (53) 

B. Differential reaction rates for neutrino emission 

As can be seen from Eqs. (|38|l . (|48|l . and H50|l . the cross sections (ail'n ^) and for the neutrino absorption 

processes v^, + n —^ + p and 9^ + p ^ e'^ + n, respectively, have the same generic form 



/ Q^ /' + 00 /• /'27r 
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where K = K or K' , and T is the integrand involving the relevant amplitude squared and distribution functions. If 
we use the differential reaction rates with respect to cos Q^, to describe the neutrino emission processes e~ +p —f Ve + 'n- 
and e'^ + n ^ i>e+p, then these rates also have the generic form in Eq. H54|l . This follows from the symmetry between 
the forward and reverse processes. In particular, the transition amplitudes squared and |7^+„| are identical to 

|7^^„|^ and \T[,^p\^, respectively. By taking advantage of the symmetry between the neutrino absorption and emission 
processes, numerical computation of the cross sections for the former and the differential reaction rates for the latter 
is greatly simplified. 

For e~ + p ^ Ve + n occurring in the material above the protoneutron star, the differential reaction rate is 
dAl-J =^ ^ /"Ldfc,. 1 1 



dcose, - J Mkp^,^P^Sp)dkp, E / 2^ L3 e(^=/^)-"e + 1 



-3^3, 1 1-7- |2 



47r2 ^ J (27r)3 L-^L-^ tL^ 

where rji, is the degeneracy parameter characterizing the Fermi-Dirac distribution function of the electrons. Integrating 
over the neutron momentum as in Eq. (|37|l . we obtain 



dA^ ' „„„2 /) ^ 1 ^-eB/m„T ,-- 

^^e-p Up COS Oc mn i ~ 6 ' " / 



dcose^ 327r4 ,j2^f^ ei^^^^ + e-^fivB/T)-(eB/m^T) 

oo oo /•+00 !■ 



Tip—O n„—0 



.27r -{E^-m^)/T 



Similarly, we obtain the differential reaction rate for e+ + n — > i/g + p as 

d^i+'f^ ^ f ^ r LdK, 1 1 



^ J /„(k„,s„)d^fc„ J2 J 



dcose, ^ J " „ ^„ 7 2n L3 e(Se/T)+r,e + i 

E 



47r2 2^ J 271 L-^L-^ tL^ 

Sp ,TLp ifp 

Glcos'ec m^eB 1 H E"" dE 

16^3 (27rm„T)3/2 eA-nS/i^ + e-A'.B/T - - 

EE/ dfce. EE/ dfc.^ 

g-(£;,.-m„)/T 

e(i^./T)+,. + i >^^^pd$„. (57) 



X 



IV. RATES OF NEUTRINO PROCESSES IN SUPERNOVAE 



We now calculate the rates of the neutrino absorption and emission processes in Eq. for the supernova en- 
vironment near a protoneutron star that possesses a strong magnetic field. A wide range of heavy elements may 
be produced during the ejection of the material above the protoneutron star. As mentioned in the introduction, a 
key parameter governing this production is the neutron-to-proton ratio of the material 0, Q , which depends on the 
competition between the neutrino absorption and emission processes 0] . We will calculate the rates of these processes 
in the context of heavy element nucleosynthesis, for which the accuracy of these rates is especially important. In 
this context, the material above the protoneutron star is characterized by temperatures of T ~ 1 MeV, entropies of 
S ^ 100 (in units of Boltzmann constant per nucleon), and electron fractions of Yg < 0.5. For these conditions, the 
nucleons in the material are nonrelativistic and nondegenerate while the are relativistic and have a small degen- 
eracy parameter < Tje 1- The thermal distribution functions of the nucleons and have been given in Sec. IIIII 
The neutrinos emitted from the protoneutron star are not in thermal equilibrium with the overlying material and 
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their energy distribution functions are taken to be of the form in Eq. (|51|l . As discussed in Ref . 131 , Pauh blocking 
for the final states of the neutrino processes above the protoneutron star is unimportant and will be ignored. 



A. Neutrino absorption 



At a radius R above the protoneutron star, the rate of neutrino absorption per nucleon can be estimated as 



N 



^ 49.7 



IQSi erg s ^ 



10 MeV 



(ft/ at) 



10 km 
R 



(58) 



where L^, is the neutrino luminosity and has a typical value of ~ 10^^ erg in the supernova epoch of interest. The 
key quantity (a^N) in the above equation is obtained by averaging ct^n over the neutrino energy spectrum. We first 
compare various approximations for a^jy as functions of the neutrino energy E^. 

The cross sections for neutrino absorption on nucleons in a magnetic field have been derived to ©(I/tojv) as crj^j^ 
in Sec. lIII To 0(1), the zeroth order in l/niN, the cross sections are .13j 



(0,S) 
vN 



P + 3.g2 



cos 6;, + 2x 



(/T.g) 



where 



(0) 



27r 



G p cos^ Oc , ^2 



n,=0 



2npeB 



27r 

= i;. ± A, 

'(El^^f 



ir + 3g')eB 



E. 



(0) 



(£;^°^)2-m2 



2eB 



In the above equations and elsewhere in this subsection, the upper sign is for z^e + n — > e 
for + p —^ + n. In Eq. (|5| 



X 



exp(^B/r) - exp(-/xB/T) 
exp(/iS/r) + exp(-^B/T) 



(59) 



(60) 

(61) 

(62) 
(63) 

p and the lower sign is 
(64) 



is the polarization of nucleon spin, where fi is the magnetic moment of the relevant nucleon with fip — 2.79 fiN and 
Hn = — 1.91/XAr. For the case of interest here, \ijl\B/T <^ 1, so 



X = ^ = 3.15 X 10-2 



B 



MAT J \10^^ G 



(^) 



(65) 



The term proportional to (Tb,2 in Eq- (|59|l arises because the ground Landau level of has only one spin state 
while any other level has two. In Eq. (|60|l for aB,i, the product of eB and the sum gives the total phase space of the 
in the final state. In the limit ne,max 3> 1, the summation of the Landau levels can be replaced by integration and 
<7b,i approaches 



G% cos2 9c 



Ei'^JiEi'^r~ml 



(66) 



which is the cross section to 0(1) in the absence of any magnetic field. In the same limit, aB.2 is negligible compared 
with <jb,i, so cr^']^^^ approaches cr^']|(l + e^), where 



2(/±g)g 



cos 81, 



(67) 
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FIG. 1; (Color online) Comparison of various approximations for a^jv [solid curves: cr^^^'; dotted curves: iJ"!^^^ \ short-dashed 
curves: (j^'^(l-|-e;^); dot-dashed curves: cr^'f^(l-l-e^); long-dashed curves: (y^^pj (X^^xA- The dot-dashed and long-dashed curves 
for (Ju^n are indistinguishable. The results for Ov^n are shown in (a)-(c) for cos 0^ = —1, 0, and 1, respectively, while those 
for (Jv^p are shown in (d) for cosSj^ = (oa^^ has little angular dependence). These results are calculated using B = 10^^ G, 
T = 2 MeV for the matter temperature, Xn = —0.03, and Xp = 0.04. 



results from the polarization of the initial nucleon spin by the magnetic field. 

For numerical examples of the cross sections, we take B = 10^^ G. The cross sections a'v^n^ for cosOj/ = —1, 0, 
and 1 as functions of Ey are shown as the dotted curves in Figs, Chi andCJ:, respectively. The angle-dependent 
terms in Eq. H59|l for of^^^ are proportional to the difference between / and g. As the numerical values of / and g 

are close, these terms are very small. So wc only show the cross section of^^^ for cos 9^ = as the dotted curve in 
Fig-Di- All the dotted curves in Fig.^havc spikes superposed on a general trend. The varying heights of these spikes 
are artifacts of the plotting tool: all the spikes should have been infinitely high as they correspond to "resonances" 
eB, for which a new Landau level opens up. These singularities are integrable and do not give 
infinite probabilities of interaction in practice. For example, at a given E^, the thermal motion of the absorbing 
nucleons will produce a range of Eg and the cross section obtained from integration over this range will be finite. 
Thus, the spikes in cr^*]^^'' will be smeared out by the thermal motion of the absorbing nucleons, which is similar to 
the Doppler broadening of the photon absorption lines in the solar light spectrum. The effects of such motion are of 
©(I/toat) and have been taken into account by the cross sections cr^^^' derived in Sec. 1111 M Using T ~ 2 MeV for 
illustration, we show as the solid curves in Fig.^ It can be seen that where spikes occur in o'^i^pf\ there are 

only smooth bumps in cr^^^''. Clearly, cr'^^^^ is more physical than (j'^^^^ ■ 

Two more aspects of Fig. ^ require discussion. First, the bumps in cr^^^' diminish as E^ increases and become 
invisible for Ey ^ ^/ eB ^ 8 MeV. This is expected from the correspondence principle: when a number of Landau 
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levels for and protons can be occupied, the quantum effects of the magnetic field are small. As noted in Sec. IIIBI 
the absorbing proton in I'^.+p + n can occupy many levels for T ~ 1 MeV. However, for the e"*" in +p —fc'^ + n 
and the and the proton in i^e + n ^ +p, occupation of many levels requires Ei, '3> VeB ^ 8 MeV (see Sees. Ill Al 
and IIIBji . Second, while the general trends of the dotted curves for crl^'rf'' appear to follow the corresponding solid 
curves for ai],'n^ , the general trend of the dotted curve for (J^^p ^ deviates substantially from the corresponding solid 
curve for a'^l'p'. This concerns the effects of weak magnetism and recoil of the final-state nucleons, both of which 
are of 0{l/mN) and are taken into account by cr'"^^^ but not by Figure shows that these effects give small 

corrections to cr^Jn ^ but much larger corrections to crp|!jf ^ . 

To better understand the effects of weak magnetism and recoil of the final-state nucleons, we make use of the 
correspondence principle. As noted above, the effects of Landau levels become negligible for E^, ^ VeB ~ 8 MeV. 
In this case, the only surviving quantum effect of the magnetic field is polarization of the initial nucleon spin, which 
gives rise to a dependence on cos for the cross sections due to parity violation of weak interaction. Thus, allowing 
for this surviving effect, we should recover the results for no magnetic field in the limit of high E,^. In the absence of 
any field, the cross sections cr^^ to ©(l/rnjv) is (see, e.g., Refs. 0,13) 

where we have ignored terms like m'^/E'^ and A/mjv- The above zero-field cross sections assume that the initial 
nucleon spin is unpolarized, and therefore, do not depend on cosG^. If a small polarization x is artificially imposed, 
the modified cross sections should have an additional factor 1 -I- e^. The term proportional to _E^/r7ijv in Eq. H68|) 
represents the effects of weak magnetism and recoil of the final-state nucleons. The coefficient in this term is 1.01 for 
I'e + n ^ e~ + p and —7.21 for De + P ^ + n. Therefore, over the range E^ ~ 10-50 MeV typical of supernova 
neutrinos, the correction due to the above effects is ~ 1-5% for the former reaction but amounts to ~ —7% to —36% 
for the latter reaction. The importance of these corrections has been discussed in other contexts p3.l25j . 

Using Xn = —0.03 and Xp = 0-04 corresponding to i? = 10^^ G and T = 2 MeV, we show cr^']|(l -|- e^^) and 
<jI^^{1 + e^) as the short-dashed and dot-dashed curves, respectively, in Fig. ^ The small increase from (t1°^ to ailli 
and the much larger decrease from <T^Jp to cr^^Jp given in Eq. (|68|l can be seen from this figure. In addition, as expected 
from the correspondence principle, the general trends of the dotted curves for ul^if •* closely follow the short-dashed 
curves for ajyf^ni^ + ^x^) ^ ^0 MeV and the solid curves for crf/^jf ^ become indistinguishable from the dot-dashed 

curves for ajy^li{l + e^^) in the same regime (see Figs.^-c). However, while the relation between the dotted curve for 
o^^p^ and the short-dashed curve for (J^^p{l -t-e^p) is in accordance with the correspondence principle, the solid curve 
for CTp^jf^ clearly stays above the dot-dashed curve for crp^p(l + e^p) at E,, > 25 MeV (see Fig. ^). This apparent 
violation of the correspondence principle for cTp^p^' and (Tp^p(l -f e^^) is caused by the slightly different treatments of 
the reaction kinematics in calculating Cp^jf ^ and (Xp^p. 

We have used the transition amplitudes to 0(1/771 at) in calculating both a^^^^ and cr^^. However, we have treated 
the reaction kinematics exactly for cr^^^^ [assuming nonrelativistic nucleons, see Eqs. H40H44|l ] but only to 0{l/mN) 

for crf,^. This difference does not affect the comparison between ail'n^ and aj^'Xi^ + Exn) ^s the total correction from 
weak magnetism and nucleon recoil is small to 0{l/mN) in this case, and the terms of orders higher than 0(l/rnjv) 
are even smaller. In contrast, the importance of the weak magnetism and recoil effects for + p ^ e'^ + n enables 
terms of orders higher than ©(l/mAr) to give rather large corrections to the cross section. Such terms are included in 
(J^^p^ due to exact treatment of the reaction kinematics but not in Up^^p. In Ref. [26j |. the zero-field cross sections for 
neutrino absorption on nucleons were derived to ©(l/m^v) but with reaction kinematics treated exactly. We denote 
these cross sections as cr^^V ■ For consistency with the rest of the paper, we ignore radiative corrections and the effect 
of the Coulomb interaction between the final-state particles for Ve + n + p, both of which were taken into 

account in Ref. [2^. It was shown in this reference that (7l^*J is more accurate than (ip^p. We show <j'"J'^\l + e^) ^s 
the long-dashed curves in Fig. ^ It can be seen that the long-dashed curves for ai 

(1 + ^Xn) are indistinguishable 

from the corresponding dot-dashed curves for (IH- e^^) over the range of E^ shown but the long-dashed curve for 
CTp^*J{l + e^p) lies significantly above the dot-dashed curve for <jl^]p{l + e^p) ^ 25 MeV. In addition, the solid 

curves for ail'n^ and dp^p^^ settle down to the corresponding long-dashed curves for ailn + Cxn) ^nd a^*p {1 + e^p) 

(1 B) (1 B) 

at El, > 20 and 25 MeV, respectively. Thus, the cross sections aljn and a^^'p calculated above are in full agreement 
with the correspondence principle. 
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TABLE I: Comparison of various approximations for {a^N) (in units of 10 cm ^). These results are calculated using 
= 11 MeV, = 16MeV, B = 10^'' G, T = 2MeV for the matter temperature, Xn = -0.03, and Xp = 0.04. 













+ p — > e"*" + n 




cos Qy 


-1 





1 


-1 





1 



(cr,^^>(l + ex) 1-67 1.62 1.57 2.49 2.48 2.46 

(cr^'^)(l + ex) 1-70 1.65 1.60 2.06 2.04 2.02 

(cr^,'^'>(l + ex) 1-69 1.65 1.60 2.11 2.09 2.07 

{a^°f^) 1.65 1.57 1.50 2.50 2.46 2.42 

{a'-^jf'^) 1.68 1.61 1.54 2.11 2.09 2.07 



We now calculate the average cross sections {a^N) using the neutrino energy spectra in Eq. if^ . We take -qi,^ — 
jjp^ = 3, {Ey^) = 11 MeV, and {E;^^) = 16 MeV. For these parameters, T^^ = 2.75 MeV and Tp^ = 4 MeV. Adopting 

the same B, T, and Xp as for Fig.[ll we give (ct^']|)(1 + e^)' + ^x)' (^iw^)!! + ^x)' ( 

for cos6p — —1, 0, and 1, respectively, in Tabled As discussed above, cr^^ and cr^^^ differ from a'f^^ due to the 
effects of weak magnetism and recoil of the final-state nucleons. These effects slightly increase the cross sections for 
i/g + n — > e~ + p but substantially decrease those for + p — > e+ + n. As can be seen from Table n (cri,il)(l + t^^) 
and {a^*n){\ + e^^) are only a few percent larger than (crl°k)(l + ^xr,) but (crp^p)(l + e^^) and (cTp^p )(1 + e^p) are 

~ 20% smaller than (CTp°p)(l + e^p)- ^he differences between (ct^tv^'') and {(j'f^'') are similar. On the other hand, 
the effects of the magnetic field on the average cross sections are small for both Vi,+n-^e~+p and + P — > e+ + n. 
For B = 10^^ G assumed above, (crf-^if ^) is at most 4% smaller than (crl^n)(l + ^Xn) or (o'^-Jn )(1 + Cxn) while {crl^'p ^) 
is indistinguishable from {crl^*J){i + ^xp)- This is because with {E^J — 11 MeV and {Ej^J = 16 MeV, the important 
energy range for determining the average cross sections has E^, > ^/eB ^ 8 MeV, for which the effects of Landau 

f 1 -B) f 1 B) 

levels are small. As {Ej^J is substantially larger than {E^J, the magnetic field affects (tXp^p ) even less than {ahjn )■ 
The results for absorption of supernova neutrinos on nucleons can be summarized as follows. Generally speaking, 
one can use (cr^^) — (cr^*]!) to estimate the corrections without magnetic fields, and use (o'^jv^^) — (o'I^n) estimate 
the corrections due to magnetic fields. For B < 10^^ G, en) 
with an accuracy of 1%. For the same field strength, the effects of magnetic fields are not important for {a^^p), 

and ( good estimate of {<y\^^p^) with an accuracy of ~ 1%. Note that we have ignored radiative corrections 

(see, e.g., p3l) and the effect of the Coulomb field of the proton on the electron wave function (see, e.g., H^). These 
factors give corrections at the level of ~ 2% 26]. To account for them, we suggest calculating (crle") and {(J^*J) as 
in Ref. and use the results in the above estimates for (Crieff ^) and (fTp^jf 



(1 s) 

is a good estimate of (aljn ) 



B. Neutrino emission 



The differential reaction rates with respect to cos 0^ for + p ^ v,, + n and e"*" + n — s- + P in a magnetic field 
have been derived to 0{l/misi) as dA^^'^ ^/dcos8,y and dA|,+^Vdcos8,y, respectively, in Sec. IIIIBl To 0(1), the 
zeroth order in l/niN, the differential reaction rates are 



dA 



(0,S) 
eN 



eB 



dcos8i, 27r2 



dr^"/Vdcose, 
exp[{Ee/T)TVe] + l 



dfcp 



(69) 



where 



dr 



(0,S) 
eN 



dcos8i, 



.(0) 
eN 

2 



l + 2x 
r(o) 

^ eN 



(/Tg).g 
P 



35= 

„2 



cos8y 



P + 35- 



■ cos Qi, 



(0) _ Gj, cos^ Oc 
27r 



o,, (/±g)g 

^/2 + 352 



(/2+3g2)(i?eTA)2 



(70) 
(71) 
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FIG. 2: (Color online) Comparison of various approximations for dAeiv/dcos 0;^ [filled circles with error bars: dA^^y /dcosOi^; 

dashed lines: dX^^j!f^^ /dcosQ^; dotted lines: dA^'^vJ^/d cos 0^ — dot-dashed lines: dA^^'/dcosQ^ = A^^'/^]- The error 

bars on the filled circles represent the accuracy of the numerical calculation and the solid lines are linear fits to the filled circles. 
These results are calculated using B = 10^^ G,T — 2 MeV for the matter temperature, and rje = 0. 



In the above equations and elsewhere in this subsection, the upper sign is iov e~ + p ^ v,, + n and the lower sign is 
for + n ^ T>e + P- In Eq. H70() . (5n^,o is the Kronecker delta. For comparison, in the absence of any magnetic field, 
the differential reaction rates to 0(1) are 

dA(°) r r(°) d^h 



dcose^ J exp[(£;e/T)=F?/e] + l (27r)3' 

We have also calculated dA^^Vdcos 9,y to 0(1/771 jv) using the prescription in Ref. Note that both dA^'jlJ/dcosGj, 
and dA^'^*Vd 

cos©;/ are independent of cosO^y. 

To compare dA^^^^VdcosG;, with dA^^*^^ '/d cos 9;,, we take B = 10^^ G, T 2 MeV, and rje ^ 0. The differential 

reaction rates dA^j^y /dcosO^, are numerically calculated for ~ 0, 7r/4, 7r/2, 37r/4, and vr and shown as the filled 
circles with error bars in Fig. |21 Here and elsewhere in this subsection, the error bars for our results represent the 

f 1 B) 

accuracy of the numerical calculation. To very good approximation, the rates dA).^ /dcosO^^ are linear functions 
of cos Qi, as shown by the solid lines fitted to the numerical results in Fig. El The rates dxfj^^"^ /d cos as func- 
tions of cosGjy are shown as the dashed lines in the same figure. It can be seen that relative to dA^^'^Vdcos 6;^, 

e p 

dA|;_'p /dcos9;y is smaller for cosG^^ > 0.15 but larger for cosO^^ < 0.15 due to corrections of 0{l/mjf). So we 

expect that when integrated over cosO^y, the difference between xf-^^ and A^.^'^'' is small. In contrast, corrections 

of C'(l/r7ijv) make dA^+'fVd cos0;y smaller than dAg+^'' /d cos for all values of cosG^^. As discussed in the case of 
neutrino absorption, such corrections are due to the effects of weak magnetism and recoil of the final-state nucleons, 
which tend to affect the processes involving De more than those involving i>e. These corrections can also be seen by 
comparing the zero-field results dA^*^ /d cos and dA^'^*Vd cosQ^, which are shown as the dotted and dot-dashed 
lines, respectively, in Fig. [3 Note that for the parameters adopted above, the magnetic field decreases the rates for 
+ p ^ ly,, + n but increases those for e'^ + n ^ 9^ + p. 

To further explore the effects of the magnetic field on the rates of neutrino emission, we consider two representative 
sets of supernova conditions: (T, 5", Ye) = (2 MeV, 50, 0.5) and (1 MeV, 100, 0.5) for cases I and II, respectively. 
For each case, the electron degeneracy parameter rj,, can be obtained from the equations of state as discussed in 

Ref. nil. We calculate the rates A^^^^ for a number of values of B (4 x lO^^^-l.e x lO^^ G for case I and 2 x lO^^-lO^^ 

G for case II) and show the results as the filled circles with error bars in Fig. |3| The corresponding rates xfj^^'^ as 

functions of B and the zero-field results A^^ and A);^ are shown as the dashed curves and dotted and dot-dashed 

lines, respectively, in the same figure. As can be seen from Fig.|31 in the limit of small B, the dashed curves for X^^j^^"^ 

agree with the dotted curves for A^'jlJ. The approach of A^ + ^'' (filled circles) to the zero-field limit A^+*' (dot-dashed 
line) is also clearly demonstrated for case I. As a large number of Landau levels must be included in the calculation 
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FIG. 3: (Color online) Comparison of various approximations for Aejv [filled circles with error bars: A^,^^'; dashed curves: 
-^iiv^'; dotted lines: A^'J^; dot-dashed lines: A^^']. These results are calculated using two sets of conditions: (T, S, Ye) = 
(2 MeV, 50, 0.5) for (a) and (b) and (1 MeV, 100, 0.5) for (c) and (d), where T is the matter temperature. The error bars 
on the filled circles represent the accuracy of the numerical calculation. 

for small i3, it becomes computationally prohibitive to demonstrate the behavior of A).^ in this limit to the fullest 
extent. Nevertheless, the relation between A^,^^^ and ^^^j^'' for small B clearly agrees with that between A^,^^ and 

a(0) 

The dependences on B for the rates of neutrino emission shown in Fig. O require discussion. The effects of the 
magnetic field on these rates have been noted for the specific case of i? = 10^^ G, T = 2 MeV, and = shown in 
Fig. 121 More generally. Fig. |31 shows that A^.'^'^'' and A^,^'^'' first decrease with increasing B to reach some minimum 

values and then increase with B. In contrast, Ag+^'' and a|,+^'' appear to increase monotonically with B. The above 
results can be understood by considering two different effects of the magnetic field on e^. On the one hand, a stronger 
field confines more to the ground Landau level, thus reducing the average energy. This tends to decrease the 
rates of neutrino emission. On the other hand, a magnetic field changes the phase space according to 



d^/ce e.B 



(27r)3 47r2 



/ + 00 
dfcez. (73) 



Thus, the phase space increases with B, which tends to increase the rates of neutrino emission due to the increase 
in the number density of e^. The competition between the above two factors then determines the dependences on B 
for the rates of neutrino emission. 

To show quantitatively the two effects of the magnetic field on discussed above, we compare the average energy 
(i?e)_B and the number density {pe)B of in a field with the corresponding quantities for no field, (Ef.) and pe, 
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FIG. 4: (Color online) The ratios l^Ee) b / {E^) (a) and (pe)B/pe (b) as functions of B (solid curves: T — 1 MeV for the matter 
temperature; dashed curves: T — 2 MeV). The limiting case where all are in the ground Landau level is shown as the 
dot-dashed {T — 1 MeV) and dotted (T = 2 MeV) lines. These results are calculated using rje = for simplicity. 



respectively, for a wide range of B in Fig. 0] As < ^ 1 for the supernova conditions represented by cases I 
and II, we take rje — for simplicity. The major difference between these two cases lies in the temperature. The 
ratios {E^) b / {Ee) as functions of B for T — 1 and 2 MeV (cases II and I) are shown as the solid and dashed curves, 
respectively, in Fig. The corresponding ratios (pe)s/Pe are shown in Fig. For large B, it is appropriate to 
consider the limiting case where all are in the ground Landau level, and therefore, {Ef.) b / {E^) is a constant and 
{Pe)Bl Pe increases linearly with B. These limits are shown as the dot-dashed and dotted lines for T = 1 and 2 
MeV, respectively, in Fig.^] As can be seen from this figure, {E^) sj {Ee) monotonically decreases with increasing B, 
eventually approaching the constant limit, while {pe)Bl pe monotonically increases with B, eventually approaching 
the limiting linear trend. The combined result of the two effects is that XeN decreases with increasing B in weak 
field regime, and starts to increase in strong field regime after some turn-over point. From dimensional analysis, 
we expect the field at the turn-over point to be B^ ~ E'^^je with ii^ofi being some typical energy of the particles 
participating in the reaction. Because of the threshold, e~ participating in e~ -I- p — > I'e + is more energetic than 
e+ in e+ -I- n ^ +p. So Be is larger for panels (a) and (c) than for panels (b) and (d), respectively, in Fig. 13 The 
turn-over points correspond to ~ 2 x 10^^ G in panel (b) and B^ < 10^^ G in panel (d). However, the turn-over 
in these two panels is much weaker than that in panels (a) and (c) so that x!f+^^ and Ag^^' appear to increase 
monotonically with B for B > 10^^ G. In addition, because E^s is higher for higher T, Be is larger for panels (a) and 
(b) (T = 2 MeV) than for panels (c) and (d) (T = 1 MeV) in Fig. 01 

In summary, we note that the rates XeN are sensitive to the temperature T of the supernova environment regardless 
of B: lowering T by a factor of two reduces X^-p by factors of ^ 30-50 and Ae+„ by factors of ~ 6-20 (see Fig. 
In contrast, the average cross sections {a^^) only have minor dependence on T (mainly through the polarization of 
nucleon spin) so that the rates X^n essentially scale with the radius R as X^n oc R~^. For the temperature profile 
in the supernova environment of interest, the rates X^n dominate X^n- Therefore, so long as the former rates are 
calculated accurately, the latter can be estimated using X^^^^ to good approximation for B < 10^^ G. 



V. CONCLUSIONS 

In a previous paper [T^ . we calculated the rates oi v^, + n ^ -\- p and i/g + p ^ e+ + ri in supernova environments 
with strong magnetic fields assuming that the nucleon mass mjv is infinite. We also applied these rates to discuss the 
implications of such fields for supernova dynamics. In the present paper, we have taken into account the effects of a 
finite niN and developed a numerical method for calculating the above rates to ©(l/m^r) for similar environments. 
Rates with such an accuracy are required for application to supernova nucleosynthesis. 

We have shown that our results have the correct behavior in the limit of high neutrino energy or small magnetic 
field. We find that for typical supernova Ve energy distributions, magnetic fields of i? '--^ 10^^ G reduce the rate of 
i/g + n — > e~ -|- p while the ©(l/m^v) corrections due to weak magnetism and nucleon recoil increase this rate. These 
two opposite effects tend to cancel. On the other hand, the reduction of the rate of i^e + P ^ e+ + n by the 0(l/mjv) 



17 



corrections dominates the magnetic field effects for B < 10^^ G and typical supernova De energy distributions. We 
also find that for typical supernova conditions relevant for heavy element nucleosynthesis, the rates of e~ +p —^i/f. + n 
and + n — > 9(, + P first decrease and then increase with increasing B. As it is extremely time consuming to 
numerically calculate to ©(I/tojv) the rates for the above processes in strong magnetic fields, we recommend that for 
B < 10^^ G, the following approximations be implemented in models of supernova nucleosynthesis. For Vg^+n —> e~ +p, 
it is simple to calculate the average cross section including the magnetic field effects but no ©(l/m^v) corrections 

[(crS^rf'') in Tabled] or vice versa [{ailli) in Table By comparing the two with (ctI^ti), one can estimate the 
effects of magnetic fields and the ©(l/m^r) corrections, respectively. With these two kinds of corrections combined, 

{o'i],n) + (cl^rf'') — (fl^Ti) agrees with the result of the full calculation at the level of ~ 1%. For i>e + P ^ e'^ + n, 
the magnetic field effects on the average cross section can be ignored but the 0{l/mN) corrections should be included 
with an exact treatment of the reaction kinematics [{o'l]*J) in Table Q]. While we have ignored radiative corrections 
and the effect of the Coulomb field of the proton on the electron wave function, these factors can be included in {al]}n) 
and (cTp^p ) following Ref. |2^]. For e~ + p ^ v^ + n and e+ + n ^ Ve+P-, the rates including the magnetic field effects 
but no ©(I/tojv) corrections [Eqs. l|^ - (|7T|l ] are sufficient. 

In conclusion, we note that the cross sections of neutrino absorption on nucleons are relevant not only for supernova 
nucleosynthesis but also for determining the thermal decoupling of and i>e from the protoneutron star. For example, 
a decrease in Up^p would enable to emerge from deeper and hotter regions of the protoneutron star, thus increasing 
the average Ve energy. Accurate neutrino energy spectra are essential to models of supernova nucleosynthesis 0. 
However, our results on the cross sections for neutrino absorption cannot be applied directly to the discussion of 
neutrino decoupling from a strongly-magnetized protoneutron star because the conditions (e.g., temperature and 
density) inside such a star are very different from those considered here. We also note that neutrino scattering on 
plays a significant role in supernova explosion 12^. Similar to the case of capture on nucleons, the rates of neutrino 
scattering on above the protoneutron star will be modified substantially by strong magnetic fields. These issues 
remain to be explored in detail by future studies. 
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APPENDIX A: SPECIAL FUNCTION 

The special function In,r{C) can be written as 

InAC) = xl^.e-^'^^'^-'-^'^L^-'-iO. (Al) 
V n\ 

where i"(a;) is the generalized Laguerre polynomial defined as [2^ 

1 d" 

Ll{x) = -e":E-" — (e-"a;"+") (A2a) 



! da; 



m=0 



^ ^ ' \ n ^ m. I 'm\ ^ ' 



n-m I ml 



To calculate /n,r(C) efficiently, we use its properties given below. 

• Mirror relation 
Based on the identity 

(_l)"-'-^-("-'^)g-«(^) ^ Q«--(C), (A3) 

where 

Qr'-(C)EEr!Lr'-(C), (A4) 
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it is straightforward to show that 

InAC) = (-!)""'■/,■,« (C). (A5) 

• Recursion relations 
Using the recursion relation of the generahzed Laguerre polynomial |29| 

ir'(c) = ^^:(c)-^^i(c), (A6) 

one can show that 

-In^lAO + J-In-l.r-liO- (A7) 

n ' y n 

Using this recursion relation and the mirror relation in Eq. I|A5(I . one can prove another recursion relation 

InAC) = -\l^ln.r-l{0 + -^/n- (C) • (AS) 

Starting from the definition 

/o,o(C) = e-^/^ (A9) 
one can use the recursion relation in Eq. (IA7|) to obtain 

/ /"n 

InAO = \/^^.0(C). (AlO) 

V n\ 

Using the above result and the mirror relation in Eq. lA5p . one has 

loAO = (-1)^/^/0,0(0- (All) 
V r! 

The function ln,r{0 with n > and r > can be calculated as follows: 

1. Compute /„_i,o(C) and /„.o(C) from Eq. IjAlOp . Set r' — 1. 

2. Compute In.r'iC) from the recursion relation in Eq. (|A8p . 

3. If r' = r, finish. Otherwise, compute /„_i^r'(C) from the recursion relation in Eq. ljA7|l . 

4. Advance r' by unity and return to step|21 

APPENDIX B: REDUCED AMPLITUDE SQUARED 

The reduced amplitude squared W^^n ioy: i^e + n ^ e^ + p is defined in Eq. The amphtude Tli,^n [Eq- ISSl] 

contained in Wi,^n can be simplified using 

Jo Jo 

.,(n.-r-e-np+r,)g-,K-r,-„,+.,)0„^^^^^^(^^/2eS)/,^,,^(u.i/2eB), (Bl) 

where 0^, is the azimuthal angle of . The above result follows from j3 1^ l2^ 

g.w^.x^-i(n-r)0 = 27ri"-'"e-''("-'')'^- J„_^(wiO (B2) 



27r 
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and 



J(n~r}-in'-r'){'2V'^0ln',r'{u)I„,r{u) du = /„,„' (O^r.r' (C) : 



where Jn(C) is the Bessel function. 
Noting that 



(B3) 



(B4) 



and using Eqs. (|ll|l and (|Bip . we are able to derive the foUowing expUcit expressions of W^^n with the help of 
Mathmatica® : 

(W..„)..=i,.„=i = {f + gf{l+Vez){l + cose,)ll^^^{wl/2eB) 

+ (/ - 5)'(1 - vez){l - C08e,)ll_,„(wl/2eB) 



+2{f - .g^)— cos sin e^/„^-i,„^ (wi/2eB)/„^,„^ (w^ 

+ — I f - (/ + 5)^(1 + Vez){l + COS e^)(/c„^ + kp,) 

niN L 



- (/ + 5) (2/ + /2)(1 + t^e^) sine,A:„, 

+ /2(/ + .g)(l + Wez) sinO^Wa: 

+ (/ + 9)iV + /2)(1 + COS e.)-^ X ll(wl/2eB) 



if - .9)^(1 - Vez){l - COSe^){knz + kp^) 
+ /2(/ - .g)(l - Vez)smQ^knx 

- if - .g)(2/ + f2)il-vez) sinQ^w, 

- /2(/ - .g)(l - cose,)-^ X ll_,(wl/2eB) 



E, 



- /2(/ + .g)(l + Vez)smQ^ COS 4>w\/'2neeB 



+ if - .9) (2/ + /2)(1 - Wez) sin 61, cos (f)^^/2neeB 
+ 2{-,f + gif + /2) + /(/ - .g + /2) cos 6,) 



X C0S($„ - 



\/2neeBkn± 



(/ + .g)(2/ + /2)(l + cose,) 



y/2n^eBw± 



\/2neeBw^ 



Ee 



+ /2(/-.g)(i-cose,) 

X In.-i.n, «/2eS)/„^,„^ iwl/2eB) } 



(B5a) 



(W.,„)s,=M„=-i = V(i + „^^)(l_cose,)/2^,„J</2ei?) 

2 
ruN 

+ 2gif + /2)(1 + Vez)il - COSe,)(fce. " k,,) 



2 |- r 

1 -5(/ + 5 + /2)(l + We2)sine^A;„ 



2ti p R 

+ gif + g + /2)(1 - cose.)-^ x (u;i/2eS) 



Ee 

+ gif -g + /2)(1 + Vez)sme,w,ll„Jwl/2eB) 



+ 



gif - g + /2)fcn± cos($„ - ~ gif + g + f2)wA_ 



X (1 - COS e,) ^-^^^/„, _i,„^ («;i/2eS)/„^,„^ iwl/2eB) 

-C/e 
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- 9{f - 9 + /2)(1 + Vez) sin 6;. cos 4>w^/^^n^ 
X /„^_i,„^_i(w;i/2eB)7„^,„^_i(w;i/2eB)}, 



(B5b) 



2 

Wat 



{ -5(/ + 5 + /2)(l-Ve2)sine^fc„ 



25(/ + /2)(1 - I'e.)!! + COSe,)(fcez - fc.,) 



+ 5(/ + 5 + /2)(l + cose,) 



2neeB 



+ -9 + /2)(1 - «e.) sine„7«,/2^_i^„^K/2ei?) 
+ 9{f - 9 + f2)kn± cos(<I>„ - (pyj) - g{f + g + f2)w± 



X (1 + cose.)^^5^^"e-i,n,-i«/2eB)7„,,„,_i«/2eB) 



" 9{f — 9 + /2)(1 — Vez) sin6i/ cos^ 
x7„^_i,„^(w;i/2eB)7„^,„^(«;i/2eB)}, 



(B5c) 



(/ + gf{l - Vez){l - cose.)72^_i,„^_i(«;i/2eB) 
+(/ - gf{l + vez)il + cose.)ll^^^_,{wl/2eB) 

+2{f-9')- 



cos </)u) sin @i/Ing — l,np — l 



+ 



Ee 

{f+9f{'^ - Vez){l - COSe^){knz + kpz) 



- (/ + g}{2f + /2)(1 - Vez) sine,fc„^ 

+ /2(/ + - ■f^ez) sine^W;:, 

+ (/ + 5)(2/ + /2)(l-cose.)- 



Ee 



X 

^ ^ne — l,np 



.,{wl/2eB) 



+ 



- (/ - 5)2(1 + Vez)il + COsQ^){knz + kp 

-if- 9){V + /2)(1 + Vez) sinewy;. 



-/2(/-ff)(i + cose,) 



2neeBi 



lLnp-i{</'^eB) 



hif + 9){l- Vez ) sin Qi, cos y/2n^ 



+ (/ - 5)(2/ + /2)(1 + Wez)sine,. cos0^A/2n^eB 
- 2(/2 - 5(/ + /2) + /(/ - 5 + /2) cos e.) 
\/2n^eBkn± 



X cos($„ 



Ee 



(/ + 5)(2/ + /2)(l-cose. 



\/2neeBwi_ 



Ee 



+ /2(/-fl)(l + cose,) 



\/2neeBw± 



] 



In the above equations, Ve 



x7„^_i,„^_i«/2eB)7„,,„^_i«/2eB)}. 
A;ez / Ee ■ The reduced amplitude squared Wp^p for Pe + P 



(B5d) 
can be obtained 
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from these equations by making the substitution given in Eq. H49|l . 
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